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Q ■ Abstract 

CO \ The diagram technique allows one to calculate the correction fac- 

tors which can be used in Monte Carlo simulation of some processes. 
This is equivalent to the calculation with accounting for all or some 
part of the interference contributions. The example is presented for 
the simplest case of inelastic deuteron-deuteron interactions. 
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We will consider the discussed problem for the concrete case of intranu- 
clear cascade model with Monte Carlo simulation of events. This model is 
rather popular until now, especially at not very high energies |jj Q . It is well- 
known that all interference contributions are lost in such simulation because 
in the Monte Carlo method we can add probabilities but not amplitudes. In 
many cases it leads to not so large errors because the interference contri- 
butions are not dominated. However until now it was not possible even to 
estimate their role qualitatively. 

In the present paper we will show that it is possible, as a minimum in some 
special cases, to calculate the correction factors. Use of them is equivalent 
to the account, as a minimum, some part of interference contributions. 

Let us consider for simplicity the deuteron-deuteron interaction at energy 
a little smaller than 1 GeV per nucleon. In this situation only one secondary 
pion can be produced in each nucleon-nucleon interaction, and a secondary 
nucleon after the first inelastic collision practically can not produce another 
pion because it has no enough energy. So the main source of pion production 
in the considered case is the process of one-nucleon pair inelastic interaction, 
Fig. la. Double-nucleon pair interaction has smaller probability, the process 
of Fig. lb gives some correction to the cross section of one pion production 
whereas the process of Fig. lc is qualitatively different because it leads to 
the two pion production in one event. 

There exist also a lot of processes with elastic rescattering of secondary 
nucleons and pions but they can not change the pion multiplicity (except of 
the case of pion absorption by secondary nucleon). 

Let us consider now the processes of Fig. 1 from the point of view of 
unitarity condition. The modulo squared amplitude of Fig. la is shown 
as a cut of elastic dd scattering amplitude in Fig. 2a. If the cross section 
determined by the imaginary part of the amplitude Fig. 2a is equal to Ai, 
the cross section of the process Fig. la is equal to 
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(we neglect the difference in pp, pn and nn cross sections for simplicity). 

The modulo squared amplitudes of Fig. lb and lc correspond to the cuts 
of another diagram of dd elastic scattering amplitude which are shown in Fig. 
2b and 2c. The only difference between them is that in the case of Fig. lb 
we should take one cutted nucleon-nucleon blob with inelastic intermediate 
state and another one with elastic iViV scattering, whereas in the case of Fig. 
lc the inelastic intermediate states in both cutted blobs should be taken. So 
if the contribution of the diagram Fig. 2b (2c) to the total dd cross section 
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is equal to A 2 , the cross sections of the processes Fig. lb and lc are 
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respectively. Factor two in both these Eqs. come from the AGK cutting rules 
[||, f|] and another factor two in Eq. (2) comes from combinatoric. 

The diagram Fig. lc can not interfere with another diagrams of Fig. 1 
because it contain two pions in the final state. However the diagrams of Fig. 
la and lb can interfere that corresponds to the intermediate state of elastic 
dd amplitude shown in Fig. 2c. In accordance with AGK cutting rules this 
contribution to cross section is 
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This cross section can be calculated inside the Monte Carlo code via the 
value of 0\ c (or the correspondent number of events) : 
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So we should multiply all distributions, histograms, multiplicities, etc., 
coming from the sum of events from Fig. la and Fig. lb processes by the 
factor 

R = aia + aih ~ aia > lb < 1 (6) 

Via + 016 

and only after that add the events from the processes of Fig. lc. In particular 
one can see that the mean multiplicity of produced pions will increase because 
we add smaller number of one-pion events with the same number of two-pion 
events. 

The similar calculations can be fulfilled with the help of the same AGK 
cutting rules for more realistic cases of hadron-nucleus and nucleus-nucleus 
interactions and possibly in some another cases. Of course there exist many 
another interference contributions connected, say, with final state interac- 
tions, etc. which will be not accounted by the similar way. However some- 
times these contributions can be not essential. So one can see that the com- 
bination of Monte Carlo code which allow one to calculate, say, some angular 
distributions of produced pions, and AGK cutting rules gives the possibility 
to increase the accuracy of calculations. Possibly the similar approach can 
be used in another cases where Monte Carlo simulations are used. 
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Figure captions 

Fig. 1. Diagrams for pion production in not high energy deuteron- 
deuteron interactions. 

Fig. 2. The intermediate states of elastic dd amplitude which correspond 
(a, b and c) to the modulo squares of the diagrams of Fig. 1 and (d) to the 
interference of amplitudes Fig. la and lb. 
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